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The dynamics of heavy quarkonium systems in the strong coupling regime reduces to a quantum 
mechanical problem with a number of potentials which may be organized in powers of 1/m, m 
being the heavy quark mass. The potentials must be calculated non-perturbatively, for instance in 
lattice QCD. It is well known that the long distance behavior of the static (l/m°) potential is well 
reproduced by an effective string theory. We show that this effective string theory, if correct, should 
also reproduce the long distance behavior of all 1/m suppressed potentials. We demonstrate the 
practical usefulness of this result by finding a suitable parameterization of the recently calculated 
1/m potential. We also calculate the 1/m 2 velocity dependent and spin dependent potentials. Once 
Poincare invariance is implemented, the shapes of most of the spin independent potentials are fully 
predicted in terms of the string tension, and the shapes of the spin dependent ones in terms of a 
single parameter. 
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Heavy quarkonium systems have played a major 
role in our understanding of QCD (see [l[ for a re- 
view). The early successes of non-relativistic poten- 
tial models in describing the gross features of the 
spectrum, can nowadays be understood as emanat- 
ing from QCD in a particular kincmatical regime. 
The heavy quarks in the heavy quarkonium rest 
frame move slowly, with a velocity u < 1, which gen- 
erates a hierarchy of physical scales m 3> mv S> mv 2 
(1/mv is the typical size of the system and mv 2 the 
typical binding energy) in addition to Aqcd, the 
typical hadronic scale. This hierarchy is most con- 
veniently exploited using the effective field theories 
(EFT) of NRQCD 2, 3] and pNRQCD 0,111, which 
are built in such a way that they are equivalent to 
QCD in the kinematical regime they hold (see Q 
for a review). It was shown in [5j that in the case 
mv ~ Aqcd the relevant degrees of freedom of pN- 
RQCD (and hence of QCD) reduce to those of non- 
relativistic potential models. The potentials to be 
input in pNRQCD, however, have precise formulas 
in terms of objects computable from QCD. Some of 
these formulas were known since long @, @, @] , but 
others were uncovered when formulating this prob- 
lem in the EFT framework, like the 1/m potential 

The potentials have been computed on the lattice 
with increasing precision llj, ll2j, |l3|, |l4|, |l5 1 . Con- 



venient and economical parameterizations of lattice 
data are necessary in order to include the potentials 
as simple functions in the Schrodinger equation. For 



the static potential the naive addition of the short 
distance one-gluon exchange potential and the long 
distance linear potential, as predicted by the effec- 
tive string theory (EST) fl6i ] , which is known as the 
Cornell potential [r7|, provides a good description 
of lattice data and has been very successful in phe- 
nomenological applications. Corrections to the long 
distance linear behavior can be calculated in a sys- 
tematic manner in the EST [3, [3] (see also [2Cj|). 
For the subleading potentials, so far the only con- 
straint which has been used for such parameteriza- 
tions, is that at short distances, the potentials must 
approach their perturbative expressions. The long 
distance behavior has traditionally been a matter of 
guess work, being quite common the use of poly- 
nomials in 1/r (lately powers of r have also been 
used). The aim of this letter is to show that the 
EST also predicts the long distance behavior of the 
1/m suppressed potentials, and hence may become 
an extremely useful tool in order to find suitable pa- 
rameterizations of lattice data. 

The static potential can be obtained from the vac- 
uum expectation value of the rectangular Wilson 
loop W(T,r) [HI. The EST hypothesis maintains 



that at long distances (tAqcd 3> 1) this expecta- 
tion value can be obtained from a string action, 



lim (0\W(T,r)\0) = Z [ T>£ l e 



(1) 



where Z is an unknown constant, and = (, l (t,z), 
I = 1,2, are the transverse components of the string, 
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which fulfill the boundary conditions £ l (t,r/2) = 
£ (t, —r/2) = 0. The string action may be written 

as 

String = ~kJ dtdz (l - ld^ l d^ (2) 

where n is the string tension. This action is cor- 
rected by higher order terms in the EST counting, 
and can be obtained as a long wavelength limit of 
the Nambu-Goto action. Equations |T|) and ([2|) give 
rise to the following prediction for the long distance 
behavior of the static potential [l8j |. 

V W{r) = Kr + n-^- (3) 
12r 

where /i is an unknown constant. This result agrees 
with lattice data for k ~ 0.21 GeV 2 

The 1/m suppressed potentials are given by ex- 
pectation values of suitable operator insertions in 
the rectangular Wilson loop (see [ljj HH for con- 
crete formulas). Since the large distance behav- 
ior of the expectation value of the Wilson loop is 
given by an EST, it is natural to expect that the 
suitable operator insertions that the 1/m potentials 
need also have a representation in the EST. In or- 
der to pin down the mapping it is convenient to ex- 
press the operator insertions in a gauge invariant 
fashion. This is achieved by introducing two spin- 
less (Grassmann) fields ip and x- "0 annihilates a 
static source in the fundamental representation at 
the point r/2 = (0,0, r/2) and x creates a static 
source in the anti-fundamental representation at the 
point —r/2, {ip^ip} = {x\ x\ — 1? t ne remaining 
fixed-time (anti-)commutators being zero. The QCD 
Lagrangian is then augmented by 

SL QCD = ^(t){id - 9 A (t,r/2))ip(t) + 

xHt)(ido-gA (t,-r/2)) X (t) (4) 

The expectation value of the rectangular Wilson 
loop W(T, r) can be rewritten as 

(0|W(T,r)|0) = (0|O(|,r)Ot(_|, r) | ) ( 5 ) 

0(t,r)=xHW(t,- r -;t, r -)m (6) 

<j)(t, r; t, r') is the straight Wilson line joining the 
points r and r' at the time t. In this formalism 
the insertions of chromoelectric and chromomagnetic 
operators (see (fTg)) below and ref. [13]) correspond 
to insertions in ([5]) of the following gauge invariant 
operators, 

- x t (i)E *(t,-£) x (t) , - x \ t )W{t,- V -)x{t) (7) 



For instance, let us denote as (E'(*, r/2) E*(t', r/2)) 
the expectation value of the insertions of two chro- 
moelectric fields at the points (t, r/2) and (t', r/2) of 
the Wilson loop (T/2 > t > t' > -T/2). We have, 

<E*(t,r/2)E l (i',r/2)) = 
{Q\0{^,v)^m\t*-)m (8) 

x^E'^MOoVf^lO) 

This way of rewriting the operator insertions in the 
Wilson loop is especially convenient for the mapping 
into the EST. In the limit T — > oo, which is taken in 
the computation of the 1/m suppressed potentials, 
the chromoelectric and chromomagnetic insertions 
reduce to correlation functions of the gauge invariant 
operators ([7|). These correlation functions can now 
be mapped into the EST as correlation functions of 
some suitable EST operators. 

Therefore, what we have to do is to find a rep- 
resentation of operators like (0 in terms of string 
variables, under the guidance of the global symme- 
tries of the system. The latter correspond to the 
Dhoo group, the symmetries of a diatomic molecule 
(changing P by CP), and time reversal. In order to 
identify the implementation of the symmetry in the 
EST, it is convenient to choose a worldsheet param- 
eterization in which evolution is described by time, 
the zeroth coordinate of the string, and the labeling 
by the z coordinate, the last coordinate of the string, 
as it has already been implemented in For the 
building blocks of {7}, we have the following trans- 
formation properties with respect to the generators 
oiD hoo (z = (0,0, z)): 

• Rotations with respect to the z-axis 

E l (t,z) -> R ij W{t,-i) 
B ,; (£,z) -> FC j W(t,z) 

m - m , x (t) -> x(t) (9) 

• Reflection with respect to the zx-plane 

E*(i,z) -> p ij W(t,z) 
B%z) -» -p ij B j (t,z) 

i>(t) -> m , x (t) - x(t) (io) 

• CP 

E*(t,z)-f (E*) T (t,-z) 
B l (i,z)^-(B l ) T (t,-z) 

m^x*(t),x(t)^r(t) (ii) 
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Under time reversal they transform as follows: 
• T 



E l (i,z) -> E l (-t,z) 
B 4 (t,z) -» -B 4 (-t,z) 
^(t) - V'(-i) , X(t) ~ 



x(-t) 



(12) 



In these equations, R l i is the rotation matrix, p % 3 = 
diag(l, — 1, 1), and T stands for transpose (with re- 
spect to color indices). On the string theory side, 
the building blocks, namely the string coordinates 
£ 4 (f, z) (with £ 3 = z), transform as follows: 



Rotations with respect to the z-axis 
Reflection with respect to the zx-plane 



• CP 



• T 



(13) 



(14) 



(15) 



(16) 



We find that the following mapping satisfies the sym- 
metry requirements, 



X Ht)B l (t,- r -) X (t) 
X t(t)E 3 (t,-I) x W 



A'e lm d t d z C(t, g) 



in (|17[) will be corrected by higher order operators 
in the EST counting. The expression of the 1/m 
potentials in the EST will be obtained by substitut- 
ing the operators on the lhs of (| 1 T[) by the operators 
on the rhs of (JXTJ) and calculating the expectation 
values with the EST action ©. 

Let us illustrate it by calculating the EST expres- 
sion of the 1 /m potential. For this potential we have 



22] 



V (1 >°\r) 
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dtt((W{t r -)W{Q V -))) c (18) 



where ((• • •)) means that the expectation value of the 
operator insertions in the Wilson loop (e.g. (8)) is 
normalized to the expectation value of the Wilson 
loop ([5]), and the subscript c stands for connected. 
Hence the EST representation is 



! A 4 



dttd z d z .G ll F {t, r ^0j~) 



(19) 

where G l j?(t, z',t, z') = (£*(t, z)£ m (t', z')). This in- 
tegral is most easily computed by performing a Wick 
rotation to imaginary time. For the calculation of 
the correlator we obtain 



G l ^{it lZ ;it',z') 



rim 



i 



in 



4lTK 

cosh [£(*-*')] +cos[^(2 



z')] 



cosh [*(*-*)] -cos[f(z-z>)] 



(20) 



The time integration in (JTSJ) suffers from an UV di- 
vergence, which may be regulated by introducing a 
cut-off for small times. The contribution from this 
cut-off is just an additive constant to the potential, 
which may be absorbed into the additive constant 
that appears in the EST result for the static poten- 
tial ([3]). Up to a constant term, we then obtain 



V 



(M)( r) = ^!ln 



(21) 



A 



(17) 



-A" 2 



Hence we obtain the non-trivial result that the 1 jm 
potential must grow logarithmically at large r. Let 
us compare this result with available lattice data. 
We fitted a curve of the form V^^ if) = alogr + b 
to the data in [3 at 3 = 6/g 2 = 6.00 and r > 0.2 
i\ t u t u z z, \l, 2' ^ * ' 2 m - Note that this range already corresponds to the 

intermediate and long distance regimes r > Aqq D . 
The result is plotted in Fig. 1. As we can see, the 
fit is very good, with a reduced chi-square x 2 l^df = 
0.93 . Of course, for phenomenological applications 



where I, m = 1,2 and A, A', A", A'" ~ Aq C d are 
unknown constants with dimension of mass. The 
assignment above agrees with the early assignment 
in ref. [23j]. The EST provides an expansion of 
the physical observables in terms of 1 /rAQCD , trans- 
verse string coordinates must be counted as 1 / Aqcd , 
whereas d z and do like 1/r. Hence the expressions 



1 We have considered the errors of the different lattice points 
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FIG. 1: The lattice data for V (1 ' 0) (r), fitted to the EST 
prediction V^ 1,0 '(r) = alogr + b. 



where c£P is a matching coefficient of the NRQCD 
Lagrangian, which is inherited by the spin-orbit in- 
teraction (see @). v£g°\r) is UV divergent and 
requires regularization and renormalization. This is 
not a problem of the EST itself but rather one inher- 
ited from the static limit of QCD. The introduction 
of the static fields ip(t) and x(t) makes the solution 
of the problem straightforward. Indeed, whenever 
we have a time ordered product of local operators, 
contact (local) terms of dimension equal or smaller 
than the sum of the dimensions of the operators 
must generically be added in order to obtain finite 
results. In the case of V^g (r), which involves the 
time ordered product 

e lm ^ (t)B% ^(f)^(0)E m (0, |)^(0) (25) 



(see for instance |25[) a short distance piece com- 
patible with perturbation theory (~ l/?~ 2 ) should be 
"added" to the long distance behavior above. 

Some of the 1/m 2 potentials are related to the 
correlator (f2T)|) , and hence can be easily obtained 
from it, 



(id), 



(22) 



The velocity dependent potentials Vjp ,0 \r) and 

V^p (r) may receive non- vanishing contributions at 
NNLO, and, hence, up to logarithmic corrections, 



they are expected to scale as (f) ~ V^' L '{r) ~ 
C/r (Vp2'°\r) may develop a constant piece due to 



a contact term, similar to the ones appearing in (|26| 
below). We obtain from (|22|) the following model- 
independent predictions for the long range behavior 
of these potentials, 



V# 0) (r) 



= -I 



r 2_d_y(l,0)( \ 

dr \ / 



= 4 (23) 



Let us next turn to the potentials involving chro- 
momagnetic fields (spin dependent potentials). We 
obtain for the spin-orbit potentials, 



T/(2,0) 
V LS 



r nr 2 



2„(1) A / A 2 



9 c> 



(24) 



only the following terms are possible 2 , 

e lm ^(t)B l (t, |)^W^( )E m (0, |)^(0) 
e lm ^(t)B L (t, ^(i)V f (0)E m (0, £ty(0) 
■ft (c 1 5'{t)+c 2 S'"(t))^(0)i>(0) 
+ ic 3 S'(t)^ (0)E 3 (0,^(0) 



(26) 



ci, C2 and C3 are real constants. The term with ci is 
subleading in the EST counting, but c\ and C3 are 
not. We use the same regularization as for y( 1 '°)(r) 
and add a suitable contact term corresponding to the 
EST representation of the terms with ci and C3 in 
(|26p , which turn out to be proportional to the iden- 
tity operator, in order to make the final expression 
finite. The coefficient /1 2 appearing in (pM)) depends 
on the finite piece of this contact term and must 
be considered an additional free parameter. For the 
spin-spin potentials we get zero at LO, which is con- 
sistent with the argument put forward in (23j . How- 
ever, at NLO they might receive non- vanishing con- 
tributions. Up to logarithmic corrections, we expect 
them to scale as V^2 ~ Vg^ ~ C/r 5 , which may 
explain the sharp drop observed in lattice calcula- 
tions [lij]. Note that these contributions would be 
to 2 /Aq CD enhanced with respect to the one found 
in [H]. 

Before closing, it is interesting to explore the con- 
straints that Poincare invariance imposes on the po- 
tentials [27| with regard to the EST results above. 



uncorrelated. When the correlations are taken into account 
X 2 /Ndf becomes larger but still of order one [2rl . 



2 Note that (0)i/>(0) is the identity operator in the subspacc 
spanned by i/>t(0), and hence operators involving higher 
powers of it are redundant. 
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The Gronies relation 
tion [29] fix £ h 



and the first BBMP rela- 



and A 2 to 



Me = «/2 , 5 A 2 



(27) 



The other BBMP relations are satisfied without any 
further constraints. This is a remarkable result. It 
fixes the coefficients of the 1/m potential and of the 
velocity dependent potentials in terms of the slope 
of the static potential (the string tension k). For 
the 1/m potential the fit value of the coefficient 
a = 0.095 GeV 2 , whereas the previous relation gives 
a = k/tt = 0.067 GeV 2 . The difference may be 
due to two reasons: (i) the lattice data of [24[ are 
not in the continuum and, hence, small violations 
of Poincare invariance are expected, and (ii) higher 
order terms in the EST, which have not been con- 
sidered, the most important of which goes like C/r 2 , 
up to logarithms. 

In summary, we have shown how EST can be used 
to extract the long distance part of the 1/m sup- 
pressed potentials. As an example, we have quan- 
titatively compared with lattice data in the case 
of the 1/m potential and have found an excellent 
agreement. We expect a similar agreement for the 
remaining potentials. When Poincare invariance is 
used, the shapes of the spin-independent potentials 
are fully predicted (at LO in the EST expansion), 
and the shapes of the spin-dependent ones are given 
in terms of a single parameter. 

We believe our results are important from two dif- 
ferent points of view. On the one hand, we have 
obtained for the first time a satisfactory parame- 
terization of the 1/m potential at long distances, 
which can now be used to compute the 1/m cor- 
rection to the heavy quarkonium spectrum. On the 
other hand, there is no available proof of the idea 
that QCD is equivalent to EST at long distances. 
Our results provide a number of new ways to test 
whether this idea is valid or not. 
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